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1 INTRODUCTION 



The X-ray transient GRS 1915+105 was discovered by the 
GRANAT satellite in 1992 (Castro- Tirado et al. 1994). 
It has since shown erratic activity with recurrent rises to 
a maximum X-ray luminosity of 3x10^^ ergs/s (Chaty et 
al. 1996, hereafter C96). VLA observations at the time 
of outburst led to the discovery of relativistic ejections of 
plasma clouds with apparent superluminal motions, possibly 
a smaller-scale analogue to the jets observed in active galac- 
tic nuclei and quasars (Mirabel & Rodriguez 1994, hereafter 
MR94). GRS 1915-1-105 has a highly variable spectral in- 
dex, a hard X-ray tail at >150 keV, and often exceeds the 
Eddington luminosity limit for a neutron star. Due to these 
characteristics and its similarity to the other Galactic su- 
perluminal source, GRO J1655-40 (with a compact object 
mass of 7Mq; Orosz & Bailyn 1997), GRS 1915-^105 is a 
suspected black hole binary. Radio observations show that 
it is at a kinematic distance D — 12.5±1.5 kpc from the 
Sun; combined with measurements of the hydrogen column 
density along the line of sight, a visual extinction of Av = 
26.5±1 mag has been determined (C96). Variable radio and 
infrared counterparts to the X-ray source have been found 
(Mirabel et al. 1994). However, the optical counterpart has 
only been detected at /, with a magnitude of 23.4; upper 
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limits on the B, V, and R magnitude are 25.9, 26.1, and 
26.1 respectively (Boer, Greiner, & Motch 1996). 

The extremely high visual extinction ensures that GRS 
1915+105 is impossible to study optically. Therefore, in an 
attempt to constrain the nature of the binary system, an 
ongoing program of infrared observations of GRS 1915+105 
is currently underway. Of particular interest is the na- 
ture of the mass-donating star. On the basis of its spec- 
tral morphology and the calculated absolute K magnitude, 
Castro-Tirado et al. (1996, hereafter CT96) suggested GRS 
1915+105 to be a low-mass X-ray binary (LMXB). However, 
further K-hand spectroscopy of the counterpart revealed 
several prominent emission features but no detectable ab- 
sorption features which would be indicative of a late-type 
secondary. On the contrary, the -ft'-band spectrum is very 
similar to those of high-mass X-ray binaries (HMXB) with a 
Be star as the mass-losing component (Mirabel et al. 1997, 
hereafter M97). In addition, M97 show that CT96 under- 
estimated the visual extinction along the line of sight to 
GRS 1915+105. Using our corrected Ay, a striking simi- 
larity becomes apparent between the observed absolute IR 
magnitudes and colours of GRS 1915+105 and those of other 
HMXBs, including SS433 (the classic stellar source of rela- 
tivistic jets; Margon 1984) and A0538-66 (Allen 1984). 

JHK photometry of the counterpart has revealed both 
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short- and long-term variability; changes of one magnitude 
over an interval of 24 hours and two magnitudes over sev- 
eral months have been observed, but no periodicity has been 
seen (C96). However, all previous IR photometry of GRS 
1915+105 has been either short term, spanning only a few 
days, or highly erratic, with a handful of observations scat- 
tered across an extended time period. In order to search 
for the long-term periodicity (~15-45 days) that would be 
expected in an eccentric HMXB, a high-resolution (~1 day) 
light curve with a baseline of several months is necessary. To 
this end, we observed GRS 1915-1-105 in the 7^-band from 
July to September 1996, producing the first long-term IR 
light curve of this source. 

Since the discovery of the superluminal ejections, GRS 
1915+105 has been under almost continuous surveillance 
at X-ray and radio wavelengths. The Rossi X-ray Timing 
Explorer (RXTE) and the Burst and Transient Source Ex- 
periment (BATSE) on the Compton Gamma Ray Observa- 
tory (CGRO) have monitored the X-ray behaviour of GRS 
1915+105 extensively (see e.g. Greiner et al. 1996, here- 
after G96, and Harmon et al. 1996), while numerous ra- 
dio observations have been performed using the Very Large 
Array (VLA), Nancay Telescope, Green Bank Interferome- 
ter (GBI), and Ryle Telescope (Rodriguez & Mirabel 1997; 
Rodriguez et al. 1995; Pooley 1996). At these wavelengths, 
GRS 1915+105 shows fluctuations of variable amplitude on 
a variety of timescales. As in the IR, no consistent periodic- 
ity has been found. However, some X-ray and radio observa- 
tions have shown a quasi-periodicity of ~30 days (Rodriguez 
& Mirabel 1995). Numerous attempts have been made to 
correlate the X-ray and radio activity of GRS 1915+105, 
but the wide range of behaviour exhibited by the source has 
made such comparison complex. 

Here we present multi-wavelength light curves of GRS 
1915+105 from the period July to September 1996, includ- 
ing ii'-band photometry from the 1.8m telescope at Lowell 
Observatory, RXTE 2-10 keV, BATSE 20-200keV, and 1414 
MHz and 3310 MHz observations from the Nancay radio 
telescope. By exploring the source's behaviour simultane- 
ously at multiple wavelengths, we attempt to place some 
constraints on the nature of the system. In addition, we 
discuss possible interpretations for the correlations (or lack 
thereof) between the light curves. This study is the first to 
include extensive IR photometry as well as X-ray and radio 
information in a multi-wavelength study of GRS 1915+105. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Infrared 

Our K-haxiA images were obtained over a total of 28 nights 
during observing runs in 1996 July, August, and Septem- 
ber using the Ohio State Infrared Imager/Spectrometer 
(OSIRIS; DePoy et al. 1993) on the Perkins 1.8m telescope 
of the Ohio State and Ohio Wesleyan Universities at Lowell 
Observatory. The images were taken with the f/7 camera 
which provides a 2'7 field of view at a resolution of 0'.'63 
pixel"^; our seeing was typically I'.'S. On each night, a series 
of five 30-second (coadded) exposures were obtained, with 
offsets of ~ 20" between each exposure. The images were 
processed by first running a simple interpolation program 



Table 1. K magnitudes of GRS 1915+105 



UT Date 


K magnitude 


# of images 


Comments 


1996 




(sets of 5) 




2/7 


12.42+0.05 


2 




3/7 


12.37+0.09 


3 




5/7 


13.09+0.30 


1 


mostly cloudy 


8/7 


13.37+0.25 


2 


mostly cloudy 


11/7 


13.11+0.04 






12/7 


13.49+0.12 




variable clouds 


13/7 


13.29+0.15 




variable clouds 


14/7 


13.43+0.08 






15/7 


13.48+0.04 






6/8 


12.95+0.06 






7/8 


12.76+0.20 




poor weather 


8/8 


13.14+0.06 






10/8 


13.00+0.09 






11/8 


13.20+0.08 






12/8 


13.14+0.05 






13/8 


13.23+0.09 






14/8 


13.19+0.10 






16/8 


13.45+0.18 




poor weather 


20/8 


13.48+0.14 






27/8 


13.37+0.09 






4/9 


13.51+0.08 






8/9 


13.72+0.05 






9/9 


13.47+0.15 






15/9 


12.85+0.06 






16/9 


13.17+0.02 






20/9 


12.93+0.38 




poor weather 


22/9 


12.99+0.18 




poor weather 


23/9 


13.00+0.02 







to remove bad pixels. A median-combined sky image was 
then created from the five exposures and subtracted from 
each image. The resulting sky-subtracted frames were then 
fiat-fielded with an image constructed from the difference of 
dome flats obtained with the flat-field lamps on and off. 

The data were analyzed using the IRAF reduction task 
"apphot" with a 4.5 pixel aperture. The aperture was cho- 
sen to maximize the target signal relative to the sky back- 
ground; it also excludes a star 4" northeast of the target. 
Relative photometry was then performed using several field 
standards. The July- August monsoon season caused recur- 
rent poor weather during the three month run, producing 
large photometric errors on affected nights, on the order of 
+0.3 magnitudes; however, more typical errors are ±0.08 
mag. On the night of 6 August, absolute photometry of the 
local standards was obtained and used to calculate the K 
magnitudes of GRS 1915+105 throughout the three-month 
light curve. A journal of the IR observations with K magni- 
tudes is presented in Table 1. 



2.2 Radio 

Radio monitoring of GRS 1915+105 was initiated with the 
Nancay radio telescope in late 1993 (Rodriguez et al. 1995, 
hereafter R95) at 21 cm and 9 cm wavelengths. It was in- 
terrupted in 1995 and resumed in May 1996 with some im- 
provements. In order to minimize any possible outside in- 
terference, the observations were carried out at 1404, 1410, 
1416 and 1424 MHz (21 cm) and 3263, 3290, 3335 and 3349 
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MHz (9 cm), all with 6.4 MHz bandwidths falling in pro- 
tected radio astronomy bands of HI and CH, save for 3290 
MHz. Each daily run consisted of averaging RA drift scans 
centered on the source. The Nancay radio telescope has half- 
power beamwidths of 4' x 22' at 21 cm and 2' x 10' at 9 
cm. The confusion level due to nearby sources is usually 20 
mjy and 5 mjy, respectively, but due to the proximity of 
the Galactic plane, it rises to 400 mJy and 100 mJy in the 
vicinity of GRS 1915-1-105. It is thus essential to observe the 
source at the same hour angles so that the antenna beam 
profile remains constant. 

The flux density values are obtained by subtracting a 
"reference" drift scan or template from each daily average 
drift scan and fitting a Gaussian at the source position. This 
template was built by averaging earlier days of observation 
when the source had a flux density <10 mJy in the same 
wavelength range, as measured by interferometers such as 
the GBI (Foster et al 1996) or the VLA. The flux den- 
sity scale was regularly checked against the calibrator 3C433 
with assumed flux densities of 12.40 and 5.60 Jy at 21 and 
9 cm respectively. Although the number of drift scans is not 
constant from day-to-day, the typical errors for the Nancay 
data set are 10 mJy (21 cm or 1414 MHz) and 20 mJy (9 
cm or 3310 MHz). 



2.3 RXTE 

The RXTE All-Sky Monitor (ASM) has been continuously 
observing bright X-ray sources in the range of ~2 to 10 
keV since 1996 February 20. The X-ray observations of GRS 
1915-1-105 presented herein are extracted from the "quick- 
look results" public archive provided on the World Wide 
Web by the ASM/RXTE team. The hght curves for both 
the individual 90-second exposures (typically 5-10 per day 
per source) and for daily average flux are available. We have 
chosen to use the one-day average counts for GRS 1915-1-105 
rather than the individual "dwells" as this resolution cor- 
responds to that of our IR and radio data. Additionally, 
GRS 1915-f 105 exhibits extreme X-ray variability over very 
short timescales (seconds- minutes; G96); the one day aver- 
age counts smooth out these fluctuations, better represent- 
ing the long timescale variations in which we are interested. 
For further details about the instrument and the methods 
used in the ASM data reduction and error calculations, see 
Levine et al (1996). 



2.4 BATSE 

The BATSE Large Area Detectors (LADs) have been used 
extensively and with great success to monitor hard- X-ray 
sources using the earth-occultation technique. The data are 
sampled at 2-second intervals (although the effective time 
resolution here is ~one 90- minute spacecraft orbit). The use- 
ful limiting sensitivity of the earth-occultation technique is 
typically ~100 mCrab. For a complete description of this 
technique and its capabilities see Harmon et al. (1992). 

We constructed a light curve spanning the interval 
covered by our IR campaign using the standard earth- 
occultation data products obtained under the auspices of a 
CGRO Guest Investigator Program designed to study X-ray 
novae outbursts. The fluxes are derived from summations of 



~ i to several days, using the weighted average of typically 
2-3 LADs. In this case, count rates were converted to pho- 
ton flux by assuming a power-law index of 2.8 and applying 
an absolute calibration. Energy channels covering approxi- 
mately 20-200 keV were included in the analysis. 



3 DISCUSSION 

The five three-month light curves of GRS 1915-f 105 we ob- 
tained are shown in Figure 1. 



3.1 The Infrared Light Curve 

Our /i'-band light curve appears in the top panel of Figure 
1. GRS 1915+105 exhibits variability over one magnitude, 
from K ^^12. 5 to ~13.5. The data show maxima at UT dates 
3 July, 7 August, and 15 September. A visual inspection of 
the curve between the latter two maxima indicates an appar- 
ent 40-day modulation. An attempt to search for a regular 
periodicity on this timescale was unsuccessful, as a 40-day 
period appears inconsistent with the ~30 day separation of 
the first two maxima. However, we note that it is possible 
that our observations began after a true peak; therefore a 
^40-day cycle cannot be ruled out. 

The shape of the curve and the long timescale for the 
variations are intriguingly similar to several known Be/X- 
ray binary (XRB) systems with eccentric orbits and long 
periods (generally ^ 20 days; van den Heuvel & Rappaport 
1987). The IR photometric characteristics of GRS 1915+105 
are especially similar to the LMC Be/XRB A0538-66 
(X0535-668). A0538-66 has Port =16. 6d and shows corre- 
lated optical and X-ray emission on this period when the 
system is in an active state (Skinner 1980). Both sources 
show ^ 1 magnitude variability at K (M97), and both reach 
Lx ~ lO"'^ erg/s at outburst peak (Mirabel et al. 1996; 
Charles et al. 1983). Also, similar to our GRS 1915+105 
light curve, A0538-66 shows large cycle to cycle variations 
in its precise peak positions, mimicking a variable period. 

The inclination of the jets in GRS 1915+105 is well con- 
strained to be i = 70±2°(MR94). However, this cannot be 
taken as a tight constraint of the system inclination; for ex- 
ample, the jets of GRO 1655-40 are tilted 15°with respect to 
the orbit (Orosz & Bailyn 1997). If GRS 1915+105 similarly 
has an orbital inclination tilted within +20°with respect to 
the jets, we would expect IR variability on the orbital period 
due to X-ray heating of the companion star and disk; how- 
ever, the amplitude of these variations is partially dependent 
on the spectral type of the mass-donating star. Substan- 
tial orbital flux variations are generally expected from high- 
i LMXBs, where X-ray heating dominates the light curve. 
In HMXBs, more moderate variability (~ 10-20% in the op- 
tical) is the norm (van Paradijs & McClintock 1995). We 
note, however, that A0538-66 shows >1 mag variability at 
K (Allen 1984). In addition, the flrst maximum in our light 
curve is -^0.3 mag brighter than the two subsequent max- 
ima, indicating possible short-term variability such as that 
seen by C96. It therefore seems likely that the IR variabil- 
ity in our 7^-band light curve does not result from a single 
cause, but from several different processes, perhaps with an 
underlying ~40-day modulation. 

In addition to orbital variability, there are a number of 
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possible sources of IR emission in GRS 1915+105 which may 
account for the observed changes in magnitude. As hsted 
in M97, these include (1) IR emission from the accretion 
disk (CT96), (2) free-free emission from an X-ray driven 
wind (van Paradijs et al. 1994), (3) time- variable Doppler- 
broadened spectral line emission from ions in the relativistic 
jets (M97), (4) thermal dust reverberation of energetic out- 
bursts (Mirabel et al. 1996), and (5) synchrotron emission 
from IR jets (Sams et al. 1996). Of these possibilities, (3), 
(4), and (5) are related to ejection events and would there- 
fore produce changes in the IR magnitude correlated with jet 
activity. (Note also that the existence of the reported IR jets 
is in question; see Eikenberry & Fazio 1997.) Unfortunately, 
during this interval there were no high-resolution radio ob- 
servations which would have provided conclusive informa- 
tion about ejection events which may have occurred during 
our observations, and discerning a clear X-ray or radio sig- 
nature of such events (discussed in the following section) is 
problematic. We also cannot rule out possible contamination 
of the IR light curve by emission from the accretion disk. In 
addition, Be/XRBs exhibit an IR excess (on the order of ~1 
mag at K) associated with free-free emission in the circum- 
stellar shell (see e.g. Norton et al. 1991). Variation would 
then appear during the Be star phase change (formation or 
dissipation of the shell). The timescale over which such a 
change occurs is not known, although it is likely to be long 
(^months to years). As such, it is unlikely to be the cause of 
the variability seen in our data. Finally, enhanced IR emis- 
sion could be produced as a result of advective accretion in 
the inner accretion disk (discussed in section 3.3). 

The fundamental question is, what model of GRS 
1915+105 can reproduce the characterisitics of our IR light 
curve? Any such model must account for the apparent in- 
dependence of the IR and X-ray/radio emission from the 
system, as it appears that the peaks in the IR light curve 
are not correlated to the peaks in the radio and X-ray data 
(although the first two IR maxima precede increased radio 
emission by ~10 days). Be stars undergo erratic outbursts of 
equatorial mass ejection; a compact object in orbit around 
the Be star may therefore appear as a transient hard X- 
ray source, recurrent on an orbital timescale during such an 
outburst. However, our light curve shows no evidence of an 
IR/X-ray correlation which we might expect to see as a sig- 
nature of additional X-ray reprocessing of material newly 
accreted from the cloud during the "close approach" of the 
two components (Apparao 1985). In and of itself, this is not 
sufficient to discount the Be/XRB hypothesis. First, an X- 
ray periodicity correlated with increases in IR emission at 
periastron could easily be lost within the large, erratic X- 
ray activity from the disk. Second, as in other Be/XRBs, 
this correlation may not always be apparent; during active 
states of A0538-66, the optical/X-ray correlation was not 
consistently observed (Howarth et al. 1984). 

In fact, there is compelling IR, X-ray, and radio evidence 
for a circumstellar "shell" of gaseous material surrounding 
GRS 1915+105 (see the detailed discussion in Mirabel & 
Rodriguez 1996); such material is a natural consequence of 
Be star mass ejection. The IR magnitude of the system in- 
creases during mass loss events, when gas is injected into the 
surrounding region (which may already contain remnants of 
previous ejections). As the compact object sweeps through 
the cloud of gas near the Be star, both the shell of gas and 



the accretion disk of the compact object are disrupted and 
heated, increasing IR emission, resulting in peaks in the light 
curve recurrent on an orbital timescale. Quasi-periodicity 
may develop as trailed material from previous passes of the 
compact object through the circumstellar shell may cause 
additional peaks in IR emission during subsequent orbits 
(Boyle & Walker 1986). A complication in this scenario is 
the time required for dissipation of the ejected stellar mate- 
rial and of the compact object's accretion disk, which may 
shrink as the accreted material from periastron is consumed 
by the black hole. The X-ray and radio behaviour of GRS 
1915+105 in our data indicate that disk activity was occur- 
ring throughout the three-month interval, i.e. the accretion 
disk was always present. The IR peaks could then corre- 
spond to intervals where the compact object moved through 
the clouds of material ejected from the Be star, which is 
most likely in an active outbursting state. 

Conversely, it is difficult to explain the observed IR vari- 
abffity with a LMXB model of the system. In LMXBs, IR 
modulation is expected on the orbital period due to X-ray 
heating; the amplitude of this variation is generally less than 
50% in the optical, and decreases to 20-30% in the IR (van 
Paradijs & McClintock 1995). Therefore not only is the vari- 
ability larger than can be explained in a standard LMXB, 
the expected orbital period (usually ^ 1 day) is substantially 
smaller than that seen in our light curve. Some LMXBs do 
show long-term variability due to disk precession, but the 
amplitude of such modulation is significantly less than the 
1 magnitude seen in our data (Priedhorsky & Holt 1987). 
We have not been able to reconcile the timescale of the IR 
variability and the amplitude of the modulation within an 
LMXB model for GRS 1915+105. Therefore, we believe that 
the Be/XRB model, which can explain both the character- 
istics of the IR light curve and its apparent independence 
from the X-ray/radio activity, to be a more likely scenario. 



3.2 Radio and X-ray Light Curves 

3.2.1 Description of the Data 

Radio and X-ray light curves for July- September 1996 ap- 
pear in the lower four panels of Figure 1. While the two radio 
curves show very similar variability, the RXTE and BATSE 
data are remarkably different, indicating strong spectral 
variations. In the initial 15 days, erratic variability is seen 
in the radio and soft X-rays (RXTE). The flux level then 
drops and becomes much less variable, and a period of rela- 
tive quiescence ensues. This is particularly noticeable in the 
RXTE data, where almost no variability is seen for approx- 
imately 20 days. During the same 20-day period, the hard 
X-rays (BATSE) continue to show low-amplitude fluctua- 
tions, slowly rising throughout the soft X-ray/radio quies- 
cence. Then, at day 310 (15 August), a flare occurs at all four 
wavelengths, peaking at day 311 and then gradually falling 
off over ~15 days in the hard X-rays and radio. The soft 
X-rays continue to show large amplitude fluctuations while 
the flux at the other wavelengths continues to decrease. On 
day 334, however, a large one-day spike in the radio flux oc- 
curs, with simultaneous smaller spikes in the X-ray flux. The 
radio and hard X-ray emission then continues to decay, al- 
though the BATSE flux subsequently increases substantially 
(Figure 2). The soft X-rays do not decrease simultaneously 
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with the other wavelengths, remaining strong and variable 
through 1 October (day 357). However, an examination of 
the complete one-year RXTE light curve of GRS 1915+105 
(Figure 3) shows that the soft X-rays decayed gradually from 
the September peak (day 339) over the subsequent 100 days, 
and have returned to a low, quiescent flux level. 

3.2.2 Comparison of the Light Curves 

Despite several radio and X-ray studies which have been per- 
formed on GRS 1915-f 105 and other galactic jet sources, the 
complex behaviour of GRS 1915-1-105 is not well understood, 
nor has a clear signature of jet ejection events been positively 
identified. Foster et a7. (1996; hereafter F96) found the hard 
X-ray (BATSE) and radio (2.25 and 8.3 GHz) emission to 
be generally correlated, and classified active radio emission 
into two states, "plateau" and "flaring". Plateau states are 
characterized by a total flux density level of ~100 mjy, with 
a duration of days to weeks. Flaring states are always pre- 
ceded by plateau states, have a very rapid rise time (<1 
day), and are usually related to a decrease in the hard X- 
rays during the radio decay. Their data showed the hard X- 
ray outbursts to be correlated with plateau radio states and 
occasionally with the occurrence of strong radio flares. They 
also detected a "plateau-flaring" state, characterized by er- 
ratic radio emission superimposed on an underlying plateau 
state caused by rapid (<1 day) flares, the flrst time such 
behaviour has been seen in GRS 1915-f 105. Applying these 
criteria to our radio (3.3 and 1.4 GHz) and BATSE data, it 
appears that the majority of our radio curve (days 270-338) 
can be described as a plateau-flaring hybrid state, similar to 
that seen by F96 in October-November 1995. As expected 
from the F96 classification scheme, throughout this period 
the hard X-rays are in outburst, with BATSE flux levels 
consistently above 0.05 ph/cm^/s. We also see a strong ra- 
dio flare coincident with a peak in the hard X-rays (days 
310-327), similar to (but stronger than) that seen by F96 
around MJD 49990, although the intense 1-day radio flare 
at day 334 has only a small corresponding increase in the 
BATSE flux. The F96 data show one flare (10 August 1995) 
interpreted as an ejection event due to its rapid onset and 
decay time of a few days, as was detected in April 1994 
(R95); this ejection event was conflrmed by VLA observa- 
tions (Mirabel et al 1996). The flare of days 310-325 in our 
data shares similar characteristics with both the April 1994 
and August 1995 ejections. 

Harmon et al. (1996; hereafter H96) have also at- 
tempted to deflne a hard X-ray signature for ejection 
events in GRS 1915-1-105 by examining BATSE and muhi- 
frequency radio data. In their examination of the behaviour 
of GRS 1915+105 during its 1993-94 outburst, H96 found 
an anti-correlation between the radio flux density and the 
hard X-ray flux. This is in sharp contrast to the positive 
correlation seen by F96, a change which H96 attribute to 
the difference in radio flux levels during the two periods 
of activity. During the 93-94 outburst, the radio flux was 
generally above ~100 mJy, with extended intervals substan- 
tially in excess of ~400 mJy. In 1995, the radio flux was at 
a more quiescent level of ~10 mJy with several episodes of 
^200 mJy flares. H96 therefore assert that a positive cor- 
relation exists between radio and hard X-ray emission when 
the source exhibits a low GHz flux, while the weaker anti- 



correlation appears only at substantially higher flux den- 
sities (>200 mJy). They suggest that this indicates more 
active jet production in 93-94 and differing regimes for the 
production of radio emission during the two epochs. H96 sur- 
mise that the short-timescale anti-correlations correspond 
to ejection events (deduced from VLA interferometer data). 
Our data show some evidence for an anti-correlation between 
the radio and hard X-rays (days 270-303); the hard X-rays 
are relatively low during the radio flaring of days 272-284, 
and subsequently increase while the radio flux drops during 
days 285-308. However, in contrast to the H96 data, there is 
a clear correlation between the hard X-ray and radio light 
curves at a high flux level during the flare of days 310-327. 
For the most part, our hard X-ray and radio data exhibit 
behaviour closer to that described by F96 than H96; as the 
radio flux of our data is generally lower than that seen by 
H96, this may support their theory that the active states 
are coupled to the average radio flux density at any given 
time, thereby producing substantially different behaviour at 
different epochs. 

On the basis of the similarity in rise, duration, and de- 
cay times of the hard X-ray/radio flare during days 310- 
327 to the ejection events of April 1994 and August 1995, 
it seems likely that the flare in our data indicates that an 
ejection event occurred at this time (August 1996). If so, 
it is interesting to note that the IR light curve has a max- 
imum 8 days prior to the X-ray/radio flare. Additionally, 
the flrst IR maximum (day 267) precedes the earlier radio 
flaring (possibly anti-correlated with the hard X-rays) by 7 
days. We also note that the peak of our hard X-ray light 
curve (~0.145 ph/cm^/s — 480 mCrab), is higher than the 
peak values reported for the April 1994 (~220 mCrab) and 
August 1995 (~350 mCrab) ejection events. In our data, it 
appears that the radio and hard X-ray peaks are coincident; 
however, there is a gap in our radio data between days 311 
and 325. Hjellming (1997) noted that hard X-ray peaks are 
followed within a few days by radio peaks in ejection events 
of the superluminal source GRO J 1655-40, and therefore the 
jet ejections must begin very close to the X-ray peaks. In the 
AprU 1994 and August 1995 ejections of GRS 1915+105, the 
hard X-ray flux peaks prior to the radio. The apparent co- 
incidence of the hard X-ray and radio peaks in our data are 
not inconsistent with this pattern, as it is certainly possible 
that the actual peak during the observed flare occurred dur- 
ing the gap in our radio coverage. Therefore the pattern of 
the appearance of a hard X-ray peak followed within a few 
days by a radio peak, with subsequent correlated decays, 
may be a signature of jet ejection in GRS 1915+105. It seems 
likely, however, that this is not the only such signature; the 
anti-correlations seen by H96 (also seen in Cyg X-3; Hjellm- 
ing 1997) may also be a hallmark of ejection events. As yet 
undefined are the mechanisms which are necessary to cause 
the changes in the hard X-ray/radio behaviour surrounding 
ejection events, and hence altering the ejection signatures of 
GRS 1915+105 at various times. What does seem likely is 
that major ejection events do not occur if X-ray (hard or 
soft) active states are not accompanied by radio emission 
levels >100 mJy (M97; G96). 

In our data, there also appears to be some correlation 
of the soft X-ray light curve with the radio emission. Soft X- 
ray flaring occurs concurrently with the radio flaring of days 
274-284; the onset of soft X-ray flaring precedes that of the 



© 0000 RAS, MNRAS 000, 000-000 



6 Bandyopadhyay et al. 



radio by ~9 days, and falls to quiescence ~6 days before the 
radio flares subside. The soft X-rays also return to outburst 
simultaneously with the major hard X-ray and radio flare of 
days 310-318, which we have labelled as an ejection event. 
In the short term (~20 days after the flare), the soft X-rays 
continue to be high and variable while the radio and hard 
X-rays decline. A strong flare coincides with the 1-day radio 
flare; however, the soft X-ray emission peak appears ~5 days 
later. As seen in Figure 3, the soft X-rays then gradually de- 
cline through the subsequent months, and finally return to 
quiescence. In contrast, the hard X-rays increase to a flux 
level similar to that of the flare. Examination of the BATSE 
and RXTE light curves from 20 February - 13 November 
1996 (days 136-401) show a distinct anti-correlation of the 
hard and soft X-rays over the majority of that interval (Fig- 
ures 2, 3); only during the suspected ejection event is similar 
behaviour seen. Because of the extremely erratic variability 
in the radio and X-ray data, however, statistical tests of 
these possible correlations do not produce useful results. 

F96 note that a soft X-ray outburst in March 1995 was 
not correlated with any detectable radio emission. In an ex- 
amination of Proportional Counter Array (PCA) data from 
RXTE, G96 surmise that the highly variable soft X-ray emis- 
sion originates from instabilities in the accretion disk rather 
than from a large, hot corona. They further assert that peri- 
ods of large-amplitude X-ray variations correlated with ra- 
dio emission may be related to jet formation. This would be 
consistent with our observations, which show such a corre- 
lation during the onset of the possible ejection event in our 
data. Origination of the soft X-ray emission in the accretion 
disk may also account for both its short-term variability and 
its continued long-term activity after the radio decay, a sce- 
nario supported by application of the advection-dominated 
accretion flow (ADAF) model to GRS 1915-1-105. 

3.3 The ADAF Model 

The most common model for accretion flows is the thin disk 
model (Shakura & Sunyaev 1973), in which heat energy re- 
leased by viscous dissipation is radiated almost immediately 
by the accreting gas. While this model has been successfully 
applied to many systems, such as accreting white dwarfs and 
soft X-ray transients (SXTs) in outburst, it has been unsuc- 
cessful at describing SXTs in quiescence, particularly the 
black hole (BH) SXTs. However, in an advection-dominated 
accretion flow (ADAF) model, the radiative efficiency of the 
accreting gas is low, so most of the viscously dissipated en- 
ergy is stored within the gas as entropy, and only a small 
fraction is radiated while the remainder is advected onto 
the compact object (Narayan et al. 1997a). If the compact 
object is a black hole, then the energy disappears through 
the event horizon. As a result, the total luminosity of the 
source will be significantly less than that expected for its 
mass transfer rate. Narayan et ad. (1996, 1997b) have suc- 
cessfully modelled the BH SXTs V404 Cyg and A0620-00 
using a combination of an ADAF and thin disk model; the 
ADAF extends from the event horizon to an intermediate 
(and variable) transition radius, and a thin disk is present 
from this point until the outer disk edge. The model has 
also been applied to Nova Mus 1991 and, most recently, the 
other galactic superluminal source, GRO J1655-40 (Esin et 
ad. 1997; Hameury et al. 1997). It therefore seems evident 



that this combination of ADAF/thin disk model could use- 
fully be applied to GRS 1915-H05. 

Belloni et al (1997) have used the ADAF model to ex- 
plain rapid, regular variations seen in PCA RXTE data dur- 
ing an active state of GRS 1915-1-105 in 1996. They surmise 
that the emission comes in two distinct varieties: "constant" 
emission from the outer accretion disk (the thin disk), fluc- 
tuating on long timescales, while "varying" emission origi- 
nates in the inner disk region (the ADAF), which repeat- 
edly empties (via infall into the black hole) and refills on 
a timescale of seconds. They also note that the accretion 
rates inferred from the X-ray spectrum are anti-correlated 
to the total source luminosity, as expected in an ADAF. 
Finally, they suggest that the outbursts of GRS 1915-f 105 
are caused not by an increase in the mass transfer rate, but 
rather by a decrease from a rate high enough to smother the 
X-ray emission. In such a high state, GRS 1915+105 would 
only appear as a bright IR source. 

We have examined our data in the ADAF/thin disk 
context. If the radio and hard X-ray emission is produced 
close to the compact object, primarily within the ADAF re- 
gion of the disk, then this emission should decay following 
an ejection event, as a significant fraction of the energy falls 
through the black hole event horizon and disappears from 
view. This is borne out in our radio light curve, which clearly 
decreases to a level lower than seen at any other time in our 
data. However, the hard X-rays decline only briefly before 
rising dramatically to the earlier flare level. Any soft X-rays 
arising from the thin disk, i.e. the "constant" component, 
could continue to be strong and variable for some time, pos- 
sibly due to X-ray reprocessing. This is consistent with the 
behaviour of the soft X-rays in our data both during and 
after the suspected ejection event. As GRS 1915-H05 is in 
an active state throughout our light curve, it seems unlikely 
that we would see the anti-correlation of the IR and X-rays 
expected from a "smothered" state. With the multitude of 
other possible sources of IR emission (listed in section 3.1), 
at best we can only add this speculation to the list as yet 
another possible cause. 



4 CONCLUSIONS 

We have presented the flrst multi- wavelength study of GRS 
1915+105 to include well-sampled IR photometry. By exam- 
ining these data, we have arrived at several conclusions. 

(i) There is evidence for a periodic IR modulation on 
the order of 30-40 days. The qualitative characteristics of 
this variability are similar to those expected for a Be/X-ray 
binary in an eccentric orbit. However, it is likely that the IR 
variation arises from a combination of causes, with short- 
term IR emission superimposed on an underlying period. 

(ii) By comparing our X-ray/radio data to previously 
observed ejection events, we believe that an ejection event 
took place during the interval 15 August to 1 September 
1996, during which time we see a large simultaneous out- 
burst in the radio and both soft and hard X-rays. It is inter- 
esting to note that during this event the IR emission is at 
a low level, indicating that IR variations resulting from jet 
ejection are minimal in this instance. However, it has been 
suggested that the ejection event reported here actually be- 
gan around day 275, when the 15 GHz radio emission was 
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observed to increase substantially from its previously quies- 
cent level (Pooley & Fender 1997). If this earlier date, which 
corresponds to one of the maxima in our IR light curve, is 
correct, then jet ejection effects may contribute to the IR 
flux during this interval. 

(iii) In general agreement with H96, we surmise that 
the pattern of a hard X-ray flare rapidly followed by a radio 
flare, with subsequent correlated decays, is a signature of 
jet ejection in GRS 1915-1-105; however, the X-ray activity 
must be accompanied by radio emission levels ^ 100 mjy 
for major ejection events. We also note that this is probably 
not the only such hallmark of jet ejection. 

(iv) The qualitative behaviour of GRS 1915-1-105 in our 
radio and X-ray data is in rough agreement with that ex- 
pected from a two-component ADAF/thin disk accretion 
model. We note, however, that as yet no detailed quanti- 
tative study of the application of this model to high-mass 
X-ray transients has been performed. For example, the rel- 
ative strengths of the various types of emission may be sub- 
stantially different in a HMXB than in low-mass SXTs. In 
addition, as GRS 1915-f 105 is in a high state throughout 
our data set, the thin disk component may dominate the X- 
ray and radio emission. It is therefore difficult to draw any 
substantial conclusions from our data about the presence of 
an ADAF in GRS 1915-1-105. 

In the near future, we will be obtaining further IR pho- 
tometry of GRS 19154-105; we plan to compile a longer K- 
band light curve in an attempt to clarify the nature of the 
30-40 day variability demonstrated here. Our primary goal 
will be to identify the spectral type of the mass-donating 
star, crucial to our understanding of GRS 1915-)-105. Com- 
bining an IR light curve of that length with simultaneous 
radio and X-ray data, as we have done here, may shed fur- 
ther light on the behaviour of this complex source. 
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Figure 1: Multi-wavelength light curves of GRS 1915+105 from July-September 1996. Note the 
simultaneous flare of days 310-327 which appears in the radio and X-ray data, but which corresponds 
to a minimum in the IR light curve. Error bars in the radio and X-ray data are smaller than the 
point symbols; errors for the IR data are shown where they are larger than the symbol size. 
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Figure 2: BATSE hard X-ray (20-200 keV) light curve of GRS 1915+105, spanning the period 31 
January - 13 November 1996 (1 Crab ~ 0.3). 
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Figure 3: RXTE ASM soft X-ray (2-10 keV) hght curve of GRS 1915+105, spanning the period 20 
February 1996 - 27 February 1997 (1 Crab ~ 75 cts). 



